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Abstract 

We compute the two-point correlation functions of general quadratic oper- 
ators in the high-temperature phase of the three-dimensional 0{N) vector 
model by using field-theoretical methods. In particular, we study the small- 
and large- momentum behavior of the corresponding scaling functions, and 
give general interpolation formulae based on a dispersive approach. More- 
over, we determine the crossover exponent (j)T associated with the traceless 
tensorial quadratic field, by computing and analyzing its six-loop perturbative 
expansion in fixed dimension. We find: (j)T = 1.184(12), (pT = 1-271(21), and 
4)T = 1.40(4) for iV = 2, 3, 5 respectively. 

PACS Numbers: 05.70.Jk, 64.60.Fr, 75.40.Cx, 61.30.-v 
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Typeset using REVT^ 



I. INTRODUCTION 



In nature many physical systems undergo phase transitions belonging to universality 
classes of the O(A^) vector models. Their universal critical properties can be determined 
theoretically by considering the 0^ Hamiltonian 



n 



(Fx 



(1) 



— * 

where is an A^-component real field. Various computational methods, supported by 
renormalization-group (RG) theory, have provided accurate determinations of several uni- 
versal quantities, see, e.g., Ref. |l| for a recent comprehensive review. Among others, we 
should mention the critical exponents, the equation of state, and the correlation functions 
of the order parameter However, for some experimental systems one is also interested 

in the behavior of correlation functions describing the critical fiuctuations of secondary, 
quadratic local fields. Due to the symmetry of the theory, there are two independent quan- 
tities that are quadratic in the fundamental field (f){x): one is the local energy density 



which is O(A^) invariant; the other one is the anisotropic second-order traceless tensor 



(2) 



Tij{x) = (j)i{x)(f)j{x) - 6ij — (f){x) ■ (j){x). 



(3) 



The crossover exponent (pT associated with the traceless tensor field Tij{x) describes 
the instability of the 0(iV)-symmetric theory against anisotropy It is thus relevant 

for the description of multicritical phenomena, for instance the critical behavior near a 
bicritical point where two critical lines with 0{N) and 0{M) symmetry meet, giving rise to 
a critical theory with enlarged 0{N + M) symmetry, see, e.g., Refs. This bicritical 

behavior has been the object of new studies quite recently, since it appears in the SO{b) 
theory of superconductivity |^ , and has been observed experimentally in organic conductors 

pT| , the correlation functions Ge{x — y) = {E{x)E{y)) and 



|10| . As discussed in Ref. 
Gt{x — y) = {Tij{x)Tij[y)) are relevant in the description of strain-strain correlations in 
certain liquids and solids, where an effective coupling between the order parameter and 
the elastic deformations occurs. Moreover, in the special case N = 2, the traceless tensor 
field Tij{x) is related to the second-harmonic order parameter in density- wave systems. 



whose critical behavior belongs to the XY universality class, see, e.g., Refs. |112| , p!3| , p!l| . 
Experimentally, such behavior is observed at the nematic-smectic A transition in liquid 
crystals [|l2| , p!4| - [T9|JTl[] . In these systems the structure factor of the secondary order parameter 
Tij has been measured using x-ray scattering techniques |]TH,r^|. The crossover exponent (pT is 



also relevant ||2^ in the description of crossover effects in diluted Ising antiferromagnets with 
ra-fold degenerate ground state I^T], for instance in some diluted magnetic semiconductors 
such as Cdi_a;Mna;Te. 

In this paper we determine the crossover exponent (pT- Such a quantity has already 



been obtained in the framework of the e-expansion to three loops |22], from the analysis of 



2 



high-temperature expansions for = 2, 3, and by means of a Monte Carlo simulation 
| 23| for N = 5. Here, we consider the alternative field-theoretical (FT) method based on 



a fixed- dimension expansion in powers of the zero-momentum quartic coupling and 



perform a six-loop calculation of (pT- For the physically interesting cases = 2, 3, 5 we 
obtain 



1.184(12) 


(iV = 


2) 


1.271(21) 


{N = 


3) 


1.40(4) 


{N = 


5) 



(4) 

We also consider the correlation functions Ge{x) and Gt{x) in the high-temperature 
phase. In the critical limit, the Fourier transform GT{q) obeys a scaling law that is analogous 
to that of the fundamental correlation function, i.e. 

GT{q^t)=A^t-^-fT{q^e). (5) 

where t = {T — T^jTc is the reduced temperature, 7t = 20^ — 2 -|- a is the tensor sus- 
ceptibility exponent and ^ is the second-moment correlation length computed from the 
two-point function of the order parameter. The same scaling behavior holds for the corre- 
lation functions ^^(g, t) of systems in the Ising universality class, with a replacing 77-, i.e. 
GEilyt) = ^e^~"/£;(q'^^^)- For N >2, however, a is negative and an additional background 
term should be taken into account. In this case, in the critical limit, we have 

GE{q,t) =Be + Gs,sing(g,t) = Be + A+t-^Eiq^e)- (6) 

The background term Be is the dominant one and the singular part vanishes at criticality. 
In this case, by using positivity (unitarity in FT language) arguments, one may also show 
that < 0, as observed in experiments. 



In this paper we extend the two-loop e-expansion computation of Refs. [|Tl|,|I^. We 
compute the universal scaling functions fE^q^i"^) and /t(9^^^) using the e expansion and the 
expansion in fixed dimension d = ?>. First, we determine the small-momentum behavior to 
four loops in the fixed- dimension expansion and to three loops in e-expansion. In particular, 
we obtain accurate estimates of the experimentally relevant ratios Xe,t = ^et/^'^> where 
C,E,T is the second- moment correlation length computed from Ge,t{x) or from its singular 
part if a is negative. For instance, for iV = 1 we find 



and for iV = 2 



Xe = 0.0140(5), (7) 



Xe = -0.0017(1), Xt = 0.041(2). (8) 



Moreover, we study the large-momentum behavior of the structure factors and construct 
interpolations valid for all momenta by using the dispersive approach applied to (0(O)0(x)) 
by Bray @. 

The paper is organized as follows. In Sec. || we report the computation of the crossover 
exponent 0t to six loops in the fixed-dimension expansion and compare our results with 
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the existing theoretical and experimental estimates (Sec. [II Dl) . In Sec. |T| we report the 
computation of the structure factors. In Sec. IIl7\| we briefly summarize the expected be- 
havior of the structure factors in the critical region and set our notations. In Sec. [Ill B| 
we explain our FT calculation, whose results are presented in Sec. [Ill C[ In Sec. P!IID[ we 
finally give approximate expressions for the structure factors by using a dispersive approach. 
Appendix |^ discusses the large-momentum behavior of the structure factors. Details of the 
perturbative calculation are reported in App. O. 



II. THE CROSSOVER EXPONENT ASSOCIATED WITH THE TENSOR 

COMPOSITE FIELD 



A. Zero-momentum scaling behavior 



The zero-momentum behavior of correlation functions involving generic local operators 
0{x), such as E{x) and Tij{x), can be obtained from the free energy in the presence of an 
external field ho coupled with 0{x). Indeed, the singular part of the free energy scales as 
i 

where h is the magnetic field, and (po is the crossover exponent. Then, by differentiating 
with respect to ho, one obtains the zero- momentum correlations and the RG relations 

= 2 - a - (j)o, 

7o = -2 + a + 20O, (10) 

where the exponents (3o and 70 describe respectively the critical (singular) behavior of the 
average {0{x)) ~ and of the susceptibility xo = Ex(^(0)^(3^))c ~ 

In this section we compute the crossover exponent (pT associated with the tensor field 
Tij{x) in the fixed- dimension FT framework, by performing a six- loop perturbative expan- 
sion. Of course, the crossover exponent associated with the energy density E{x) is trivial, 
i.e. (j)E = ^ and 7^; = a. 



B. The fixed-dimension expansion: generalities 



In the fixed-dimension FT approach, one renormalizes the theory by introducing a set 
of zero-momentum conditions for the two-point and four-point one-particle irreducible cor- 
relation functions 

rSf(p) = 5.,Vh'+p' + o(p')], (11) 

rgL(O) = g I {6,,6ki + 6^^6,1 + 6uS,k) , (12) 

where e = 4 — d and d is the space dimension. They relate the mass m and the zero- 
momentum renormalized coupling g to the corresponding Hamiltonian parameters r and 
u: 
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u = nf gZu{g) Z^{g) ^ 
In addition, one introduces the function Zt that is defined by the relation 

T[]''\0) = 6,,Z,{g)-\ 



(13) 



(14) 



where T^^''^\p) is the one-particle irreducible two-point function with an insertion of 

The critical theory is obtained by setting g = g*, where g* is the nontrivial zero of the 
/3-function 



The standard critical exponents are then obtained by evaluating the RG functions 

dlnZ, 



(15) 



d\nm 
din Zt 



dlnm 



at the fixed point g*, i.e. 



Vtig) ■- 

- = 2 + r]t{g*)-r]4g*). 

V 



(16) 



(17) 



In three dimensions these RG functions are known to six loops for generic values of p6| , |27 
For = 0, 1, 2, 3, seven-loop series for ri^ and were computed in Ref. 



In order to evaluate the crossover exponent 0t associated with the operator Tjj(x), 
we define the renormalization function Z^i^g) from the one-particle irreducible two-point 
function Y^^iji) with an insertion of the operator Tjj, i.e. we set 



V%\%,,^^, = Z7,\g)A 



where 



A 



ijkl 



ilfjk 



N 



so that Zt{0) = 1. Then, we compute the RG function 



vAg) 



dlnZr 



d In Zrp 



dlnm 

and rjT = ^rid*)- Finally, the RG scaling relation 

(f)T = {2 + TjT - r]) u 

allows us to determine 



dg 



(18) 



(19) 



(20) 



(21) 
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C. The fixed-dimension expansion: six-loop results 



(1') 

We computed (0) to six loops. The calculation is rather cumbersome, since it requires 
the evaluation of 563 Feynman diagrams. We handled it with a symbolic manipulation 
program, which generates the diagrams and computes the symmetry and group factors of 



each of them. We used the numerical results compiled in Ref. p9| for the integrals associated 
with each diagram. We obtained 



VT{g) = -g- 



2(6 + A^) _3 18.312844 + 3.433275A^ - 0.21674589A^2 



-9 



, 140.79937 + 37.573408iV + 1.0362736A^2 ^ o.094342565A^3 
+9 



(22) 



-9 



(8 + iV)4 

5 1340.075 + 416.71657A^ + 17.622623iV2 - 0.91128056iV3 - 0.050833747A^^ 

(8 + iV)5 

15651.266 + 5665.6519iV + 433.68712iV2 + 1.0675503A^3 ^ o.67910559A^^ + 0.031393004A^5 



+0(9' 



NY 



where, as usual, we have introduced the rescaled coupling g defined by 

4877 



8 + 



(23) 



Field-theoretical perturbative expansions are divergent, and thus, in order to obtain accurate 



results, an appropriate resummation is required. We use the method of Ref. |Q that takes 
into account the large-order behavior of the perturbative expansion, see, e.g., Ref. |^l|. Mean 



values and error bars are computed using the algorithm of Ref. |32]. 

Given the expansion of f]T{g), we determine the perturbative expansion of (prig), Prig), 
and 7t(^), using the relations (|I0|) and (pT]). For iV = 2 we obtain |^ 0t = 1.176(4), 



1.178(3), (3t = 0.821(6), 0.825(5), and 7r = 0.355(2), 0.358(3), where, for each exponent, 
we report the estimate obtained from the direct analysis and from the analysis of the series 
of the inverse, i.e. from l/(j)T{g), etc. The two estimates obtained for each exponent agree 
within error bars, but, with the quoted errors, the scaling relations ( pUj ) are not well satisfied. 
For instance, using u = 0.67155(27) (Ref. |4[) and (3t = 0.823(6) we obtain 0r = 1.192(6), 
while using the same value of u and 7t = 0.3565(30) we have (pT = 1.1855(15). These two 
estimates are slightly higher than those obtained from the analysis of (prig) and 1/0^(5') ■ 
Clearly, the errors are somewhat underestimated, a phenomenon that is probably connected 
with the nonanalyticity p5|- p^] of the RG functions at the fixed point g*. 

In order to obtain a conservative estimate, we have thus decided to take as estimate of 
(pT the weighted average of the direct estimates and of the estimates obtained using Pt and 



7r together with the scaling relations |3^. The (very conservative) error is such to include 
all estimates. The other exponents are dealt with analogously. The final results for several 
values of N are reported in Table |. 
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TABLE I. Critical exponents associated with the tensor field Tij{x). 



N (jiT Pt 7t 

1 1.184(12) 0.830(12) 0.354(25) 

3 1.271(21) 0.863(21) 0.41(4) 

4 1.35(4) 0.90(4) 0.45(8) 

5 1.40(4) 0.90(4) 0.50(8) 
8 1.55(4) 0.94(4) 0.61(8) 
16 1.75(6) 0.98(6) 0.77(12) 



D. Comparison with previous results 

The exponent (pT can also be computed in the e expansion. Three-loop series were derived 



in Ref. 22 



iV ^^,iV3 + 24iV^ (24) 



'2(A^ + 8) 4(A^ + 



3 + 43JV4 ^ 7gg^3 ^ 3472Ar2 + 5024A^ - 48A^(5A^ + 22) (A^ + 8)C(3) 4 

8(Ar + 8)5 ^■ 

The coefficients of this series decrease rapidly; for instance, we have 

(j)T{N = 2) = 1 + O.le + 0.066^ - 0.00735899e3 + O(e^), (25) 

cP^(N = 3) = 1 + 0.136364e + 0.083959e2 + 0.0009916^ + O(e^), (26) 

for N = 2 and 3 respectively. Thus, any resummation gives estimates that do not differ 
significantly from those obtained by simply setting e = 1. For = 2, 3 we obtain 0^ ~ 1.15, 
07- ^ 1.22, in reasonable agreement — keeping into account that these are three-loop results — 
with the estimates of Table |. They are also in agreement with the estimate of Ref. 
that reports (pT = 1.16(7) for N = 2, which has been obtained by analyzing the same O(e^) 
series and the two-loop series calculated in the framework of the fixed-dimension expansion. 
The exponent 0t has been also computed in the 1/A^ expansion P^] for d = 3: 



^r-^~^o(±). (27) 

For = 16 it gives 0t = 1.80, which agrees with the FT result of Table |. 

The exponent 0t has been estimated by high-temperature expansion techniques in 
Ref. H, obtaining 0t = 1.175(15) for = 2 and 0t = 1.250(15) for A^ = 3, in agree- 
ment with the FT estimates. For A^ = 5, the exponent 0t has also been determined by 
means of a Monte Carlo simulation |2^: 0^ = 1.387(30). 



Experimental estimates of 0t are reported in Ref. We mention the experimental 

result 0r = 1.17(2) for the (2^1 + 1) bicritical point in GdAlOg [|l|]. The (3 -> 2 + 1) bi- 
critical behavior has been studied in MnF2 [0, obtaining (prp = 1.279(31). The experimental 
results obtained for a nematic-smectic-A transition reported in Ref. |18| are /3t = 0.76(4) 
and 7t = 0.41(9). 
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III. THE STRUCTURE FACTOR OF THE BILINEAR FIELDS IN THE 

HIGH-TEMPERATURE PHASE 



A. Scaling behavior 

The two-point correlation function of the fundamental field, i.e. G{x) = (0(0) ■ (p{x)), is 
of central importance because its Fourier transform G{q) is directly related to the scattering 
intensity in scattering experiments. For t — 0^, its asymptotic behavior is given by [^,0 

Giq) = C+t-y{q'e), (28) 

where C"*" is the amplitude of the magnetic susceptibility and the function f{y) is universal. 
Taking the second-moment correlation length 



^2 _ 1 kPG'(x) _ dG{q) 



2d E.G{x) ' ' dq^ 



(29) 



as length scale, the small-momentum behavior of f[y) is f{y) = 1/(1 + ?/) + 0{y'^), with very 
small 0{y^) corrections. Theoretical results for the correlation function G{q) are reviewed, 
e.g., in Ref. §. 

In this section we study the scaling behavior of the two-point correlation functions of the 
bilinear fields E{x) and Tij{x). Like the specific heat, which is given by the zero-momentum 
component of the two-point function GE{q,t), the asymptotic behavior of GE,T{q,t) for 
t — > 0"*" is not as simple as that of the fundamental two-point function. Indeed, in the 
scaling limit t — > O"*", — >■ with g^^^ fixed, RG theory predicts 

GeAq^ t) = Be,t [1 + 0(t)] + ^t-^-.-/^,T(gT) [1 + 0(t^)] , (30) 

where Be,t and A~^j, are nonuniversal constants, fE,T{y) is a universal function satisfying 
fE,T{0) = 1, and A is the exponent related to the leading irrelevant operator. As amply 



discussed in textbooks — see, e.g., Ref. |^ — the presence of the background term Be in the 
asymptotic behavior of GE{q,t) is related to the need of an additive renormalization. One 
may easily see that the same argument applies to the two-point function GT{q,t) of Tij. 

Since 77- > for all > 2, the leading behavior of the tensor two-point function is 
determined by the singular term depending on the scaling function /t(9^^^): 

GTiq,t) = A^t-^-Mq'e) [1 + O(t^) + 0{t^n] ■ (31) 

The background term gives subleading corrections of order , that turn out to be more 
relevant than the standard scaling corrections of order t^. Indeed, for the physically relevant 
cases = 2, 3, one finds that 7^ < A (A ?s 0.53 for = 2 and A ^ 0.55 for = 3, see, 
e.g., the results reviewed in Ref. [jl[). The difference decreases as A^ ^ 00, since both 7^ 
and A converge to 1 with the same 1/A^ correction. 

The same thing holds for the energy two-point function in the case of the Ising univer- 
sality class for which a is positive, a = 0.1199(7) (Ref. [^), i.e. 
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t) = Ap-'^Mq'e) [1 + O(t^) + 0{n] , (32) 



where A ^ 0.53, see e.g. Ref. On the other hand, for the 0{N) vector models with 

N > 2, since a < 0, the background term Be gives the leading behavior of the energy 
two-point function GE{q,t): 

GE{q, t) = BE + A+t-'^fEiq'e) [1 + (t^)] +0{t). (33) 

In these cases, the singular part vanishes for t = and is usually responsible for a cusp-like 
finite maximum in the specific heat at the critical point, as it is observed in experiments 
and in lattice models. This requires the nonuniversal constant to be negative (see the 
discussion in Sec. [Ill C 

In order to single out the singular behavior, one may consider the derivative with respect 
to the reduced temperature t 

WeAq, t) ^ = -^E,TA%^t-'-'^''^-WE,T{Q'^') [1 + 0{t^, f^-'-n] , (34) 

where 

WE,T{y) = IeAv) + ^yf'EAv) = 1 + 0{y) (35) 

lE,T 

is another universal function. 

1. Small-momentum behavior 
At small momentum, i.e. for y = g^^^ ^ 1, the scaling functions fE^riu) behave as 

f^{y) = l + Y.eny\ (36) 

n=l 

M^/) = l + ^a„l/^ (37) 

n=l 

Using Eq. (pS]), these expansions can be related to those of the scaling functions WE,T{y), 

WE{y) = l + J2^ny'', (38) 

n=l 

wr(l/) = l + $^anl/". (39) 

n=l 

Indeed, it is immediate to obtain 

en = en(^l + , (40) 

a„ = a„(^l + . (41) 
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Simple arguments based on perturbation theory suggest that the convergence radius Rc of 
the small-momentum expansions is determined by the two-particle cut. The singularity in 
the complex plane closest to the origin is expected to be Us = —45"^, where S^j = ^^Z^, 



2 



;ap 



and ^gap is the exponential correlation length that determines the large- distance exponential 
behavior of the fundamental two-point function. Therefore, Rc = 45*^. For the 0{N) vector 
models, is very close to one, so that Rc ~ 4. For example, S'^j = 0.999634(4) for the Ising 
universality class 5^ = 0.999592(6) for the XY universality class g, = 0.99959(4) 
for the Heisenberg universality class and S^j = l-0.004590/iV+O(l/A^2) in the large- 
limit As a consequence, for n — > cxo, 

Cn+l CL-n+l 0,n+l 1 / 
^ — ^ _ ^ ——— ^ . (42) 

4 

The constants ei and ai are related to the universal ratios Xe,t = tI ^ introduced in Refs. 
1^,11], where ^^e^t are the second- moment correlation lengths associated with the singular 



part of the energy and of the tensor two-point functions respectively. More precisely, if 
1t,e > 0, the correlation length is defined by Eq. (|29|), replacing G{q) with Ge,t{(i)- If the 
exponent is negative, then 



ei = -(G,(o)-s,)- ^^-(^^ 



The universal ratios Xe and Xt are given by Xe = — Ci and Xt = —ai. 



(43) 

g2=0 



2. Large-momentum behavior 
The large-momentum behavior of the fundamental correlation function is given by the 



Fisher-Langer formula [48 



f(y) ~ ^1 ( 1 + ^ + ^3 \ .^s 

One may derive a similar expression for the correlation functions of the bilinear fields. The 
large-momentum behavior of the structure factors can be studied by performing a short- 
distance expansion of the two-point functions Ge{x) and Gt{x). Following the method 
outlined in Ref. we obtain the corresponding asymptotic expansions for y — > oo: 

MV) « E,y-"i^^> (l + + -|2^) . (45) 

My) « T,,-/-) (l + + ^) ^ (46) 

The derivation of these formulae is reported in App. ^ Notice that for the 0{N) vector 
models with N >2, since a < 0, fE{y) increases as y — oo. 
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B. Field-theory calculations: generalities 



Because of the presence of the background term, the FT calculation of the scaling func- 
tions fE{y) and friy) requires some care. First, we define the dimensionless functions 

GeAs, y) = uGeA^, t, m), (47) 

where g is the four-point renormalized coupling. Then, in order to eliminate the constant 
additive renormalization term, we consider the derivative with respect to m of GE,T{g,y)'- 



d 

yVE,T{g,y) = m—gE,T{g,y) 



afe|(M) _ 2,«££(M) , (48) 

dg dy 



At the fixed point g*, the functions yVE,T{g, y) differ from We,t{<1, t), defined in Eq. (p4D, by 
a multiplicative factor independent of q. Therefore, the scaling functions WE,T{g,y), defined 
in Eq. (pS]), are given by 

Note that the zero-momentum functions yVE,T{g, 0) are related to the exponents '-yE,T by the 
relation 

-^=lim /?(,) ^^"^;-^(^'°^ (50) 
u g-^g* dg 



C. Field-theoretical results 

1. Small-momentum expansion 

We compute the small-momentum expansion of the structure factors to four loops in the 
fixed-dimension approach and to three loops in the e expansion. 

In the fixed-dimension approach, we first determine the expansion in powers of g of the 
coefficients Cj and di defined in Eqs. (|38| ) and (|39|). The explicit expressions are reported 
in App. |B[ In order to obtain numerical estimates we use the same resummation procedure 
outlined in the previous section. Our numerical results are presented in Tables |T| and 
PI . Note that, as expected, the ratios Ci+i/ej and di+i/di quickly approach —1/4. The 
corresponding coefficients and are obtained by using the relations ( PUD and (^l]). For 



the exponent v we use the same values reported before ||38[, while for 7-^ we use the results 
of Table |. In the case of Oj a large part of the uncertainty is due to the error in the exponent 
7t that enters the relation between a, and Oj. The results are reported in Table We 
also performed direct analyses of the coefficients Cj, Oj, considering the g-seiies that can 
be obtained from Eqs. ( ^OD and (^). The results are substantially consistent with those 
obtained by first estimating Cj and cij. In the case of they turn out to be more precise; 
we show also them in Table |I^ (third column of results). 
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TABLE II. Estimates of the coefficients Cj for several values of N. 



N 


(1 


(2/(1 






f"o/f'i 


1 


-0.170(5) 


-0.206(1) 


-0.221(1) 


-0.229(1) 


-0.234(1) 


2 


-0.155(5) 


-0.199(1) 


-0.216(2) 


-0.226(2) 


-0.232(2) 


3 


-0.142(5) 


-0.193(2) 


-0.213(2) 


-0.222(2) 


-0.230(3) 


4 


-0.133(6) 


-0.189(3) 


-0.211(2) 


-0.222(3) 


-0.228(3) 


5 


-0.126(6) 


-0.186(3) 


-0.209(3) 


-0.221(3) 


-0.228(4) 


8 


-0.111(5) 


-0.180(3) 


-0.206(3) 


-0.219(4) 


-0.227(4) 




m A "OT TTl TTT 

lABLE 111. 


Estimates of the coefficients a, for several values of N. 




N 


ai 


02/01 


a3/a'2 




05/04 


2 


-0.203(2) 


-0.224(2) 


-0.232(1) 


-0.236(1) 


-0.239(1) 


3 


-0.208(2) 


-0.226(1) 


-0.234(1) 


-0.238(1) 


-0.240(1) 


4 


-0.213(2) 


-0.228(1) 


-0.235(1) 


-0.239(1) 


-0.241(1) 


5 


-0.216(1) 


-0.230(1) 


-0.236(1) 


-0.240(1) 


-0.242(1) 


8 


-0.224(1) 


-0.235(1) 


-0.239(1) 


-0.242(1) 


-0.244(1) 


TABLE IV. Results of the coefficients Cj 


and ai for several values of N and 


from various 


analyses: (a) (d = 3) by using the fixed-dimension results for 


and Oj, and by directly analyzing 


the series for af, (b) (e-exp) by resummation of the three-loop e-expansion. 




N i 


€i {d — 3) from ei 


Oi {d = 3) from cii 


ai {d — 3) direct 


ei (e-exp) 


ai (e-exp) 


1 1 


-0.0137(2) 






-0.0145(11) 




2 


! 0.147(3)xl0-2 






0.17(2)xl0-2 






; -0.219(4)xl0-3 






-0.26(3) X 10-3 




4 


[ 0.38{l)xl0-4 






0.47(6) X 10-4 






1 -0.71(l)xl0-s 






-0.9(l)xl0-5 




2 1 


0.0017(1) 


-0.042(2) 


-0.042(1) 


0.000(2) 


-0.041(2) 


2 


! -0.017(l)xl0"^ 


0.530(3)xl0-2 


0.52(l)xl0-^ 


0.00(3) X 10-^ 


0.48(3) X 10-2 




; 0.024(2) X 10^3 


-0.85(l)xl0-3 


-0.84(3) X 10-3 


0.00(5) X 10-3 


-0.75(6) X 10-3 


4 


[ -0.041(2) X 10"* 


1.5(2) X 10-4 


1.5(l)xl0-4 


0.0(1) X 10-4 


1.3(l)xl0-4 


r 


1 0.076(4) X 10-5 


-3.0(2) X 10-5 


-3.0(2) X 10-5 


0.0(2) X 10-5 


-2.6(2) X 10-5 


3 1 


0.015(2) 


-0.047(4) 


-0.0465(8) 


0.012(2) 


-0.045(1) 


2 


! -0.14(2)xl0-2 


0.59(5) X 10-2 


0.59(l)xl0-2 


-0.13(3)xl0-2 


0.54(3) X 10-2 




; 0.19(2) X 10-3 


-0.96(9) X 10-3 


-0.96(2) X 10-3 


0.19(5) X 10-3 


-0.85(5) X 10-3 




\ -0.32(4) X 10-4 


1.8(2)xl0-4 


1.75(7) X 10-4 


-0.3(1) X 10-4 


1.5(l)xl0-4 


r 


1 0.6(1) X 10-5 


-3.4(3)xl0-5 


-3.4(2) X 10-5 


0.6(2) X 10-5 


-2.9(2) X 10-5 


4 1 


0.026(1) 


-0.05(1) 


-0.0500(6) 


0.022(2) 


-0.049(1) 


2 


! -0.23(l)xl0-2 


0.63(9)xl0-2 


0.645(5)xl0-2 


-0.22(3) X 10-2 


0.60(2)xl0-2 




; 0.31(l)xl0-3 


-1.0(2)xl0-3 


-1.06(2)xl0-3 


0.33(4)xl0-3 


-0.94(4) X 10-3 


4 


[ -0.50(2) xlO-"* 


1.9(3)xl0-4 


1.96(6) X 10-4 


-0.6(l)xlO-4 


1.7(l)xlO-4 


r 


> 0.91(3) x 10-5 


-3.7(6)xl0-5 


-3.9(2) X 10-5 


l.l(2)xl0-5 


-3.3(l)xl0-5 


5 1 


0.030(1) 


-0.053(6) 


-0.0533(4) 


0.029(2) 


-0.0528(3) 


2 


! -0.25(l)xl0-2 


0.7(l)xl0-2 


0.699(5)xl0-2 


-0.30(2)xl0-2 


0.65(2)xl0-2 


2 


; 0.34(l)xlO-3 


-1.2(2)xl0-3 


-1.16(2)xl0-3 


0.44(4) X 10-3 


-1.03(3)xl0-3 


4 


[ -0.55(2) xlO-'' 


2.1(3)xl0-4 


2.15(5)xl0-4 


-0.8(l)xl0-4 


1.85(5) X 10-4 


r 


; 1.00(3) X 10-5 


-4.2(6) X 10-5 


-4.3(2)xl0-5 


1.4(2) X 10-5 


-3.6(1) X 10-5 


8 1 


0.047(1) 


-0.060(6) 


-0.0602(1) 


0.046(2) 


-0.061(1) 


2 


! -0.35(l)xl0-2 


0.8(l)xl0-2 


0.817(5) X 10-2 


-0.43(l)xl0-2 


0.77(1) X 10-2 


2 


; 0.45(1) X 10-3 


-1.4(2) X 10-3 


-1.40(2)xl0-3 


0.62(2) X 10-3 


-1.24(2)xl0-3 




\ -0.72(2) X 10-4 


2.6(3) X 10-4 


2.58(4) X 10-4 


-1.05(3) X 10-4 


2.24(4) X 10-4 


r 


; 1.28(4) X 10-5 


-5.1(6)xl0-5 


-5.1(2)xl0-5 


1.95(5) X 10-5 


-4.4(1) X 10-5 
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TABLE V. Final estimates of the coefficients ei and ai. 



N ei ai 

1 -0.0140(5) 

2 0.0017(1) -0.041(2) 

3 0.014(3) -0.046(1) 

4 0.024(3) -0.051(2) 

5 0.030(1) -0.0533(5) 
8 0.047(1) -0.062(2) 



In e expansion we directly resum the expansions of Cj and reported in App. 0. The re- 
sults are also reported in Table |^ and are in substantial agreement with the fixed-dimension 
results. For ai, we also perform a constrained analysis that makes use of the available re- 



sults for tti in two and one dimensions. Such a method was introduced in Ref. ||5^ and 



generalized in Refs. |35,51|. In many instances it has provided quite accurate results for 



critical quantities. We use the estimates of oi in two dimensions reported in Refs. |]52| , |53 
ai = -0.0812(5), -0.1014(6), and -0.1313(9) for = 3,4,8 respectively. We also make 
use of the one-dimensional result |5^, oi = — (A^ — 1)^/(4A^^). From the analysis con- 
strained in one dimension, in three dimensions we obtain: ai = —0.0397(2) for N = 2 and 
Oi = —0.0533(3) for A^ = 5, while from the analysis constrained in two dimensions we obtain 
ai = -0.0460(3) for A^ = 3, ai = -0.0507(5) for A^ = 4, ai = -0.0621(2) for A^ = 8. Con- 
straining the analysis both in two and one dimension, we obtain: ai = —0.0458(1) for A^ = 3, 
ai = —0.0514(6) for A^ = 4, ai = —0.0625(2) for A^ = 8. These results are compatible with 
those of Table 0. 

Taking into account the above results for ei and ai, we consider as our final estimates the 
numbers reported in Table 0, where we have been rather conservative in giving the errors, 
which include all the results we have obtained. 



As already noted in Ref. [0, ei = -0/(67) + 0{e ) and ai = — 7t/(67) -|- 0(e ). These 
relations are not satisfied to order e^, see App. B. Nonetheless, they still provide very good 
approximations to Ci and Oi. For instance (see Ref. [0] for the estimates of the critical 
exponents): -0/(67) = -0.01377(8), 0.00185(10), 0.0159(8) respectively for A^ = 1,2,3, 
where the error is related to the uncertainty on the estimates of a and 7. 

The coefficient ei has also been computed for A^ = 1 by Monte Carlo simulations ||55|| . 
The numerical data are well described by — aefr(^)/(67efr(^)), where aes(t) and 7eff(^) are 
effective exponents determined from the specific heat and the susceptibility. 

It is interesting to note that the signs of and are strictly related to the signs of the 
amplitudes A~^j, and of the exponents a and jt- First, we observe that in the critical limit 
the correlation functions are nonnegative, i.e. Ge,t{x) > 0. Indeed, the lattice (p"^ model with 
nearest-neighbor couplings is exactly reflection positive and therefore, the above-reported 
inequalities are rigorously true for any value of the couplings. At criticality, they should 
hold for any model in the same universality class. Therefore, all moments are positive, i.e. 

\x\'^^Ge,t{x) > 0. If the correlation functions have the scaling forms (^) and (p2|), this 
implies 

A+T>0, (-l)"a„>0. (51) 
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For N >2, using Eq. (|33|), we obtain 

Be > 0, i-lTA^en > 0. (52) 



Relations (0) are satisfied by our results, while Eq. ( ^2]) and our result ei > imply 
< 0. Thus, although A~^rp is nonuniversal, the positivity (unitarity in FT language) of 

the theory fixes its sign. 

As a final remark, note that ai and ei are very small. The unexpectedly small value of 

ai was crucial to provide an explanation p!9|JTT[] , consistently with RG theory, of the exper- 



imental results of Ref. for density-wave systems, at the nematic-smectic-A transition 
of liquid crystals. The values of ai are quite smaller than what would be naively expected. 
The nearest singularity in the complex y plane corresponds to the two-particle cut, thus at 
large distance Gt{x) ~ Ixl"^ exp(— |x|/^y) with = ^gap/2, where ^gap is the exponential 
correlation length that determines the large-distance exponential behavior of the fundamen- 
tal two-point function. Positivity of Gt{x) would then require that the second- moment 
correlation length be smaller than As a consequence, since ^gap — ^ (see Sec. [Ill A 1| ), 
ai < 1/4. But this bound turns out to be much larger than the actual value of ai. 



2. Large-momentum expansion 

We also compute the constants Ei and Tj of the large-momentum behavior of fE,T{y)- 
Matching the large-momentum expansion of the two-loop expression of GE,T{(l,t) with 
Eqs. (|45| ) and (^6]) , we obtain 

A — N 

E, = l + ^^e + 0{e'), (53) 
i^2 = -2 + ^f^e + E,,e^ + 0(a, 



and 



Moreover, 



^3 = 2- -g^f + E3,€' + 0(6 



_ _2(4 + iV) 2(4 + iV)(20 -13iV-iV-) 

2 A-N (4-iV)2(8 + iV) ^ 23 ^ V 

2(4 + A0 (4 + iV)(56 - 34iV-iV^) , ^ , ^,^3, 



(A^=^ - 14iV2 - 140A^ - 432) - 4 2 

E22 + E32 - 2{N + 8y ^ 12(iVT8) ' 

(iV + 4)(Ar^ + 14Ar + 108) N + 4 , 
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The constants Ei and Ti are in agreement with the results of Ref. [TT[]. The divergence of 
the coefficients T2 and T3 for — 4 is related to the vanishing of the 0(e) term in the 
expansion of a 



56 



The large size of the coefficients makes it difficult to resum the perturbative series. For the 
physically interesting case of = 2 we report the result obtained by setting e = 1 and give 
as error the size of the last coefficient. In this way we obtain: Ei = 1.3(3), E2 = —0.7(1.3), 
and E3 = 0.3(1.7) for iV = 1, and Ei = 1.2(2), E2 = -0.6(1.4), E3 = 0.4(1.6), Ti = 1.6(6), 
T2 = -9(3), T3 = 8.4(2.4) for = 2. Moreover, E2 + E3 = 0.0(2) and T2 + T3 = -0.7(1) 
for N = 2. 



D. Interpolations of the structure factors 



In Ref. 



TT[] the authors discuss several approximate forms for fE,Tiy)- 

approximant for 



generalizations of the Fisher-Burford 



They present 
These approximations are 
quite crude and do not reproduce the full Fisher-Langer behavior for large y. A better 
approach based on dispersion theory was put forward by Bray [pH] . Here, we will apply the 



same method to the universal functions fE{y) and friy)- 

A generalization of the arguments presented in Ref. gives the following representation 
for friy) 



My) = 1 



sm 



7r \ 2u J 



dx Et(x) 

45+ x + y 



(56) 



where Ft{x) is the spectral function satisfying ^7^(00) = 1. We assume here that the only 
singularities of friy) in the complex plane are branch cuts on the negative real axis and 
that the leading one corresponds to the two-particle state, so that the disk \y\ < 43^ is free 
of singularities. Under this assumption, the representation (^) is exact. 

For generic Ft{x), Eq. (|56D does not give the correct Fisher-Langer behavior (^6|). 
Indeed, for |/ — cxd we obtain friy) ~ constant -|- Tiy~"''^/^'^'^\ We must thus require the 
constant to be zero. This gives the sum rule 




■00 

dxx-^-^^^^^"^ Ft{x) = 1, 

•m 



(57) 



which allows the determination of Ti once Ft{x) is given. 

Eq. ( 1^ applies also to /^(y) with the obvious replacements. However, the sum rule 
(^) requires a > and can thus be used only in the Ising case. For a < 0, Eq. (|57|) is 
replaced by 



El (Txa 
— sm — 

TT \2V 



•00 

+ / dxx~ 
'45+ 



l-a/(2i/) 



{Fe{x) 



(58) 



In order to obtain approximate expressions for the structure factors, we must assume a 



specific form for the spectral function. For this purpose, we assume, as in Ref. [^], that 
Ft{x) gives the exact Fisher-Langer behavior on the cut. Explicitly, we consider 
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1000 



FIG. 1. Universal function fE{y) obtained using Eqs. (pq ) and i 
We also report the large-y behavior, fsiy) ~ LlOOy"^'*^^"^^^ and the 
fE{y) ~ 1 - 0.013662/ + 0.001467y2 - 0.000219y3. 



for N = I. 
small-y behavior, 



Ft{x) = 1 + T2X-(^-°)/(2'^) 



7r(l — a) vrfl — a) tt'Jt 

cos z h Sm COt 



2u 2u 



IT IT -K^T 

COS h sm — cot 

2v 2v 2v 



2v 



(59) 



To completely determine the spectral function, we must specify the constants T2 and T3. 

These estimates are not very precise, 



We use here the e-expansion results of Sec. 
but the interpolation is quite insensitive on T2 and T3 separately. Indeed, what really 
matters is their sum T2 + T3 that is more accurately determined. In order to test these 
interpolations, we can compare the estimates of T\ and ai — and, analogously, of E\ and 
Cj — with those of the preceding sections. For N = 2, using T2 = —9 and T3 = 8.4, we obtain 
Ti 1.56, fli ~ —0.055, 02 ~ 0.008, which are reasonably close to the estimates reported 
before. Analogously, using E2 = —0.6 and E^, = 0.4, we obtain Ei ^ 1.00, ei ~ 0.005 and 
62 ~ —0.0007, again in reasonable agreement with previous results. In particular, the fact 
that |ei| <^ |ai| is correctly predicted by the approximation. For = 1, using E2 = —2/3, 
E3 = 1/3, we obtain Ei ^ 1.20, ei ~ —0.016, 62 ~ 0.0019, in reasonable agreement with 
what reported above. 

In Fig. |1| we report fE{y) ior N = 1 and in Figs. ^ and |] a graph of and of frill) 

for N = 2. It is interesting to note that for N = 2 the function fE{y) varies slowly and 
differs from one only for quite large values of y. Taking also into account that the pref actor 
vanishes as t ^ 0, the dependence of GE{q,t) should be hardly visible in experiments 
and in numerical Monte Carlo simulations. Moreover, in this case fsiy) > 1 for all y, so 
that, because of the inequalitites (|5^), there is an attenuation of the singular behavior for 
increasing q, as generally expected. 
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FIG. 2. Universal function fEiu) obtained using Eqs. ( pq ) and (|59[), for N = 2. 
We also report the large-y behavior, fE{v) ^ l.OOy''''^^'^^'^'^ and the small-y behavior, 
fE{y) w 1 + O.OOlTly - 0.000169?/2 + 0.0000243y3. 
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FIG. 3. Universal function friy) obtained using Eqs. (56) and (pg|), for N = 2. 
We also report the large-y behavior, /r(y) ~ 1.559?/~^'^^^^^^ and the small-y behavior, 
friy) « 1 - 0.0397y + O.OOSSy^ - 0.000852y3. 
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APPENDIX A: LARGE-MOMENTUM BEHAVIOR FOR THE BILINEAR 

CORRELATION FUNCTIONS 

In this appendix we compute the large-momentum behavior of the correlation function. 



We follow closely the discussion of Refs. |57,49| for the correlation function of the field 



1. The energy correlation function 



The basic ingredient of the calculation is the short- distance expansion of the product of 
operators E{x + y/2)E{x — y/2). For y ^ 0, see, e.g., Ref. [0, this product is equal to the 
sum of all the operators that are allowed by symmetries, multiplied by C-number coefficients, 
that take in account the short-distance behavior. The most singular contribution comes from 
the operators of smallest dimension. In this case, neglecting the contribution related to the 
identity operator, it implies 



E{x + y/2)E{x — y/2) = C{y)E{x) + less singular contributions. 



(Al) 



Now, let us consider the connected correlation function of / composite operators E{x 



. . . and its renormalized counterpart G^l^ {pi, . . . ,pi) = ZlZ^^G^^\pi, . . . ,pi). 



Then, Eq. ( [AID implies for p ^ m 



Gf (p, -p, 0, . . . , 0) ^ G{p- m)Gl-'^ (0, . . . , 0) 



(A2) 



where we have explicitly written the mass dependence of the short- distance coefficient. Since 
renormalized correlation functions scale canonically, i.e. 



Gf{p, -p,0,...,0) = m'^-2^/(p/m), 

we have 

G{p] m) = m^'^C{p/m) 
Renormalized correlation functions satisfy the Callan-Symanzik equation 

G^llipu ...,pi)= m'a{g)Gl+'\0,p,, . . . 



om og 



(A3) 



(A4) 



(A5) 



where a{g) is a RG function satisfying (7{g*) = 2 — rj, and ^72(5') = Vt{g) ~ V<i>{9)- Applying 
the Callan-Symanzik equation to the relation (|A2D we obtain, setting g = g*, 



d 1 ~ 

"^^-^2(fi'*) C(p;m)=0, 

and therefore, using Eq. (|A4|) , we have 

C{p;m) ~ m~^(p/m)~^~''^ = m~'^{p/m)~^^'^ . 



(A6) 



(A7) 
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Now, using the above-reported results and Zt/Z^ ~ m^* ~ m}/'^ ^ for g = g*, see Eq. (|16|), 
we obtain 

^G^'\p, -p) = G^'\0,0,p,-p) m'~'/^G^^\0,0,p, -p) 

^ m^~^/''C{p] m)Gf (0, 0, 0) ~ {p/m)-^'''m'^-^/'' ~ (A8) 

Integrating this equation twice with respect to t, we have 

-p) = a{p) + + ct^-^p-^/'' + o(t^-"), (A9) 

where a{p) and are unknown functions of p. Comparing this result with the scaling 
equations ( p2D and (p3|), we obtain finally Eq. (^S]). 

2. The tensor correlation function 

The calculation is analogous. The short-distance expansion of the product Tij{x)Tij{y) 
is given by 

Tij{x + y/2)Tij{x — y/2) = CT{y)E{x) + less singular contributions. (AlO) 

Now, we consider the connected correlation function with / fields E{x) and two fields 
Tij{x) with the indices summed over, Gip\pi,p2;qi, ■ ■ ■ ,qi), and its renormalized counter- 
part G^^]j^{pi,p2] qi...,qi) = Z'^Z\Z^^~'^G^^ {pi,p2\ qi ■ ■ ■ , qi). For p > m we have 

G^iP^ -p; 0, . . . , 0) ^ Ct{p; m)Gi^'^ (0, . . . , 0), (All) 

The coefficient Grip] m) scales as in Eq. ( [A4| ) and satisfies the RG equation 

CTip;m) = 0, (A12) 

where v'2 ~ Vt — V4>- Therefore 

Ct{p; m) ~ m-2(p/m) -2-2^2 +'?2 = m-\p/m)-^^+^^-''^/'' . (A13) 

As in the energy case, we consider the second derivative of G^\p, —p) with respect to t. For 
p ^ m we have 

^GPip, -p) = G^ip, -p- 0, 0) ~ m«-2/-2<^-/-4%(p, -p- 0, 0) 

^ m^-2/^-2<^^/^C(p; m)Gg^(0, 0, 0) ~ ^-(i+^t-")/^'^-!-"^ (A14) 

where we have used the fact that, for g = g*, Z^/Z^ ~ ^»?t-'?0 ^ m'^'^^'^~'^, see Eqs. (|16]), 
(pop. Integrating this equation twice with respect to t and using the scaling equation (pID, 
we obtain the large- momentum behavior (|i6|). 



•9 

m- 2?72 + V2 

om 
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(10) (11) 

FIG. 4. Feynman diagrams contributing in the four-loop computation of the two-point func- 
tions Ge,t- The black blobs indicate the insertion of the bilinear operators. 

APPENDIX B: PERTURBATIVE EXPANSION OF THE TWO-POINT 

FUNCTIONS Ge,t 

In order to compute the structure factor of the bilinear fields, we determine the one- 
particle-irriducible diagrams with insertions of two operators E or Tj^ and zero external 
legs. We use the susceptibility x as inverse mass square, so that tadpole diagrams can be 
neglected. Also, subdiagrams that correspond to diagrams of the two-point function {4>4>) 
are subtracted at zero momentum. The diagrams contributing up to four loops are drawn 
in Fig. |. The structure factors of the bilinear fields can be expanded as 

G{g,y) = uGEAn.q) = Y.{-ir'u'x"' Sfif''^ I,{q\) (Bl) 

where the sum is over the graphs without tadpoles, / is the number of loops of the graph, 
Sj the graph symmetry factor, Cf''^ the group factor, and Ij{q^) the loop integral with unit 
mass. In Table |VJ we report Sj, C^'^ . 

We computed the coefficients and di to four loops in the fixed- dimensions expansion. 
The expansion of the loop integrals Ijiq^) is reported in Table [V|. In the calculation we 
used the results of Refs. [^,^. We also used the expression of the bare coupling m as a 
function of the renormalized coupling g, 
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TABLE VI. For each diagram j contributing to the energy and tensor two-point function we 
report: the number of loops I, the symmetry factor Sj, the group factors Cj ' , and the expansion 
of the integral Ij (y) in fixed dimension d = 3. 



J 



6 
7 



I Si 



4 
4 
4 
4 



cl 

N 



_JCJ_ 
N(N-l) 



(87r)'/,-(2/) 



1 1 2 

2 2 1 

3 3 i 

4 3 I 



5 3 1 



1 

2+jV 
3 

9 

2+N 
3 



2+N 
3 

(2+jV)3 

27 , 
(2+Af)2 

9 ^ 
(2+Af)^ 

(2+AfK8+Ar) 
27 



10 4 2 (H±iXli8±AQ 



11 4 1 (2+iV)(8+jV) 
-^-^ ^ 2 27 



4 
3 
8 
9 

2(2+Af) 
3 

2(6+ Af) 
9 



16 
27 

4(2+N) 
9 

4(6+jV) 
27 

2(2+jV)(8+Ar) 
27 

8(4+ jV) 
27 



2{16+6N+N'^) 
27 



^ arctan ^ 

-0.0376821 + 0.00160802y + 0.000209975?/2 _ 0.00012236^^ 
+0.0000404108y'' - 0.00001 16689^^ + 0{y^) 

0.5 - 0.105903?/ + 0.0222193y2 - 0.00477681y^ 
+0.00104957J/'' - 0.000234591y^ + 0{y^) 

If 

hh 
hh 

-0.0266277 + 0.0012789?/ + 0.000157474?/2 _ 0.000095736?/3 

+0.000031929y^ - 9.26022 lO-^y^ + 0(/) 
0.25 - 0.0601852y + 0.0132661^^ - 0.00292294^^ 

+0.000651627y^ - 0.000147057y5 + 0{y^) 
0.322467 - 0.0786798y + 0.0175558y2 - 0.0039039^^ 

+0.0008763?/^ - 0.000198793y5 + 0(/) 
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u = m ■ 



487ig 
8 + N) 



27Ar2 + 350Ar + 1348 ^ 

1 + ^ + 5 f 

27(8 + Nf 

+ 17.3632A^2 ^ 120. 783iV + 315.831 , 



(B2) 



and the relation between x and m: 



27{N + 8y 



r -0.106993 



N + 2 



f +0{g 



Writing 



j=0 
j=0 



(B3) 

(B4) 
(B5) 



we computed the coefficients e^j- and up to j = 3. They are reported in Table |Vli . 

We computed the coefficients Cj and in e expansion to three loops, i.e. to order e^. The 
expansion of the integrals l4{y) and I^iy) was obtained by using the algebraic algorithm of 
Ref. [|g. 

We write 



N-A 



(iV + 2)(A^-4) 
(AT + 8)3 



X3X + 



X2e^ + 
{N + 2) 



N + 4 N + A 

{N + 2){N + 4) 



{N 
2 , 



X4C(3) + 



(iV + 8)5 



X5 



(B6) 



1 



(B7) 



{N + 8f (iV + 8)4 ' {N + 8f 

where A = 1.171953619344729445. The coefficients Xj are reported in Tables |V111| and ffX 



Note that lTl[| to order e^, ei ~ — 0/(67) a-^^d ai ~ — 7t/(67). These relations do not hold 
at order e^. 
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TABLE VII. Coefficients e^j- and aij, cf. Eqs. (|B|) and (| 



« 3 



1 
1 
2 
3 



eij{8 + Ny/{2 + N)^^-^^o) 



_i 

4 

1 

6 

0.0290678 
0.2941 - 0.0105012iV 



_i 

4 

1 

3 

0.0581358 - 0.0195433 iV 
0.588201 - 0.158869 - 0.044705 iV^ 



2 
1 
2 
3 



j_ 

16 
_ J_ 
15 

0.0280008 + 0.0208333 AT 
-0.046127 + 0.0223387 A^ 



16 
__2_ 

15 

0.0560016 + 0.00585959 A^ 
-0.0922537 + 0.0828459 A^ + 0.0173899 A^^ 



3 
1 
2 
3 



_ j_ 

64 
71 
3360 

-0.0179798 - 0.0114583A^ 



l_ 

64 
71 
1680 

-0.0359598 - 0.00137167A^ 



0.00517395 - 0.00217669 A^ + 0.00231481A^2 0.0103479 - 0.0333737 A^ - 0.0056229 A^^ 



1 

256 
31 
5040 

0.00733893 + 0.0044494A^ 



1 

256 

31 
2520 

0.0146779 + 0.000299779A^ 



-0.000842908 - 0.00212265Ar + 0.00162037Ar2 -0.00168582 + 0.0117321Ar + 0.00170808Ar2 



]_ 

1024 
3043 



]_ 

1024 
3043 



1774080 887040 

-0.00253604 - 0.00149678A^ -0.00507208 - 0.0000641563Ar 

-0.000390898 - 0.0015352A^ - 0.000747354Ar2 0.000781794 - 0.0037964Af - 0.000501 577 A^^ 



TABLE VIII. Expansion coefficients xi, X2 and X3 for Cj and Oj. 



Xi 



X2 



X3 



{2+N) (44-H3 Af) 
12 

-40-270 iV-42 N^+N^ 
360 

-352+1176 jV-l-180 N^-dN'^ 
10080 

1808-3072 Af-477 Af^+13 
151200 

-3392+4208 jV4-668 jV^-17 
1108800 



ei 

62 
63 
64 
65 



J_ 

12 
1_ 

120 

1 
840 
]_ 

5040 

1 

27720 



12 

-20+44 N+3N^ 
720 

-88-128 Af-9Af^ 

10080 
904+778 Af+55 Af^ 
302400 
-2544-1768 iV- 125 jV^ 
3326400 



_J_ 

72 

7 
360 
]_ 

240 

1 

1080 
]_ 

4752 



X 

72 
7 
360 
1 
240 

]_ 

1080 
1 

4752 



ai 

0-3 

05 



_ j_ 
12 

1 

120 
i_ 

840 
1 

5040 
1 

27720 
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TABLE IX. Expansion coefficients and for ej and Cj. 



ei X4 
X5 



62 X4 

X5 



X5 

64 X4 

X5 

65 X4 
X5 



-1184-348 iV-7Ar2 



32 

-{2+N) (-4112-2596 Ar-466 N^+N^) 
24 

4512+1420 N+35 
1280 

{2+N) (-170624-78048 Ar-8148 Ar2+509 N^) 
5760 

„ „ -2311264-726900 Af-17885 

'^S -^4 4587520 

586167808+604936960 Af+160060792 Af^ +8862854 Af3-1077019 Ar''+2048 

61931520 
2370848+719452 Ar+16023 
27525120 

-1006341632-1133004544 Ar-290766168 iV^-13609246 jV3+2081479 jV^-3072 

619315200 
-2525643296-724481980 Af-13348335 



155021475840 

3084090980864+3778414828800 Ar+961368439624 +41380542442 AT^ -6987518685 jV"t+2097152 

10463949619200 

~ I -1376+764 Ar+301Ar^+14 

"1 X4 

„ 7552-2784 jV- 5692 - 1667 Af ^ - 100 N* 

■^5 96 

- (3104+12300 Af+3185 Af^+UO AT^) 
02 X4 5120 

-1087744-484736 Af+26236 Af^ +32387 N'-^+2316 N"^ 
■^'5 23040 
_ „ 1445632+1966400 A^+440020 Af^+17885 

"3 X4 4587520 

568281088+262815488 AT- 1452256 Af^-13472924 Af3-1114341 Af^-9216 

61931520 

- (33881888+34540972 Ar+7018809 Af2+256368 ) 
04 X4 440401920 

-17032726784- 7719736448 Ar+90484812 AT^ +409316471 AT^ +36197264 Af'^+491520 
•^5 9909043200 

5244438496+4565788340 Af +853347495 AT^ +26696670 
"5 ■''4 310042951680 

34179580035328+14923763131008 Af-513765484268 Ar^-895119755607 AT^ -81461592926 Af'^- 1376256000 
■^5 104639496192000 
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